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A random matrix model to describe the coupling of m-fold symmetry is con- 
structed. The particular threefold case is used to analyze data on eigenfrequencies 
of elastomechanical vibration of an anisotropic quartz block. It is suggested that 
such experimental/theoretical study may supply a powerful means to discern intrin- 
sic symmetry of physical systems. 

The standard ensembles of Random Matrix Theory (RMT) [l] have had wide appli- 
cation in the description of the statistical properties of eigenvalues and eigenfunctions of 
complex many-body systems. Other ensembles have also been introduced [2|, in order 
to cover situations that depart from universality classes of RMT. One such class of en- 
sembles is the so-called Deformed Gaussian Orthogonal Ensemble (DGOE) [1, 0, 0, @] 
that proved to be particularly useful when one wants to study the breaking of a discrete 
symmetry in a many-body system such as the atomic nucleus. 

In fact, the use of spectral statistics as a probe of symmetries in physical systems 
has been a subject of intensive experimental and theoretical investigation following the 
pioneering work of Bohigas, Giannoni and Schmit [7| which showed that the quantal 
behaviour of classically chaotic systems exhibits the predictions supplied by the RMT. 
Examples of symmetry breaking in physical systems that have been studied include nuclei 
0,9, atoms 0, EH and mesoscopic devices such as quantum dots [l2| . 



In the case of nuclei, the Mitchell group at the Triangle Universities Nuclear Labo- 
ratory [1, [sl , studied the effect of isospin symmetry breaking, in odd-odd nuclei such as 
26 Al. They detected the breakdown of this important symmetry by the applications of 
two statistics: the short-range, nearest neighbor level spacing distribution (NND) and 
the long range Dyson's A-statistics [1,0]. These results were well described by a DGOE 
in which a pair of diagonal blocks is coupled. The strength of the coupling needed to 
account for the symmetry breaking can be traced to the average matrix element of the 
Coulomb interaction responsible for this discrete symmetry breaking 0, 13]. The justi- 



fication for the use of block matrices to describe the statistics of a superposition of R 
spectra with different values of the conserved quatum number can be traced to Refs. 
[l|, [14|. In the case of non-interacting spectra, i.e. if the quantum number is exactly 
conserved, the answer is a superposition of the R spectra. Since the level repulsion is 
present in each one of the R spectra, their superposition does not show this feature. 
Thus, we can say that for each spectra of states of a given value of the quantum number, 
one attaches a random matrix (GOE). For R spectra each of which has a given value of 
the conserved quantum number, one would have an R x R block diagonal matrix. Each 
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block matrix will have a dimension dictated by the number state of that spectra. If the 
quantum number is not conserved then the R x R block matrix acquires non-diagonal 
matrices that measure the degree of the breaking of the associated symmetry. This idea 
was employed by Guhr and Weidenmiiller [l3| and Hussein and Pato [3] to discuss isopin 
violation in the nucleus 26 At In reference [3J, the random block matrix model was called 
the Deformed Gaussian Orthogonal Ensemble (DGOE). 

In order to study transitions amongst universal classes of ensembles such as order- 
chaos (Poisson— >GOE), symmetry violation transitions (2GOE— >1G0E), experiments 
on physical systems are more complicated due to the difficulty of tuning the interaction 
(except, e.g. in highly excited atoms where the application of a magnetic field allows 
the study of GOE-GUE transitions). To simulate the microscopic physical systems, one 
relies on analog computers such as microwave cavities, pioneered by A. Richter and 
collaborators 15] and acoustic resonators of Ellegaard and collaborators [3, [13, 18]. It 



is worth mentioning at this point that the first to draw attention to the applicability of 
RMT to accoustic waves in physical system was Weaver (ig| . 



In the experiment of Ellegaard et al. 17] what was measured were eigenfrequencies 



of the elastomechanical vibrations of an anisotropic crystal block with a D3 point-group 
symmetry. The rectangular crystal block employed by Ellegard was so prepared as to 
have only a two-fold flip symmetry retained. Then, to all effects, the quartz specimen 
resembles a system of two three-dimensional Sinai billiards. The statistical treatment 
of the eigenfrequencies of such a block would follow that of the superposition of two 
uncoupled GOE's. 

Then, by removing octants of progressively larger radius from a corner of the crystal 
block this remnant two-fold symmetry was gradually broken. The spectral statistics show 
a transition towards fully a chaotic system as the octant radius increases. What was then 
seen was that the measured NND is compatible with a two block DGOE description but 
the A-statistics was discrepant. This discrepancy was attributed to pseudo integrable 
behavior and this explanation was later implemented with the result that the long-range 
behavior was fitted at the cost, however, of loosing the previous agreement shown by the 
NND0- 

Here we reanalyse this experiment following the simpler idea of extending the DGOE 
matrix model [E| to consider the coupling of three instead of two GOE's 0] . We show that 
within this extension both, the short- and the long-range statistics, are reasonably fitted 
suggesting that the assumption of the reduction of the complex symmetries of anisotropic 
quartz block may not be correct. Our findings have the potential of supplying very precise 
means of testing details of symmetry breaking in pysical systems. 

To define the ensembles of random matrices we are going to work with, we recall 
the construction based on the Maximum Entropy Principle [3|, that leads to a random 
Hamiltonian which can be cast into the form 

H = H + \Hi, (1) 

where the block diagonal Ho is a matrix made of m uncoupled GOE blocks and A (0 < 
A < 1) is the parameter that controls the coupling among the blocks represented by the 
Hi off-diagonal blocks. For A = 1, the Hi part completes the H part and H = H GOE . 

These two matrices H and Hi are better expressed introducing the following m pro- 
jection operators 

Pi = ^2\j><j\, (2) 

jeii 



where Jj defines the domain of variation of the row and column indexes associated with 
ith diagonal block of size Mj. Since we are specifically interested in the transition from 
a set of m uncoupled GOE's to a single GOE, we use the above projectors to generalize 
our previous model 0, 0] by writing 

rn 

H = J2 m GOE P t (3) 
1=1 

and 

m 

Hi^P^Qi (4) 

i=l 

where Qi = 1 — Pj. It is easily verified that H = H for A = 1. 

The joint probability distribution of matrix elements can be put in the form 0, 21 



P(H, a, 0) = Z N X exp (-atrH 2 - ptrHf) (5) 
with the parameter A being given in terms of a and (3 by 

\ = {l + (3/a)- 1/2 . (6) 

Statistical measures of the completely uncoupled m blocks have been derived. They 
show that level repulsion disappears which can be understood since eigenvalues from 
different blocks behave independently. In fact, as m increases the Poisson statistics are 
gradually approached. In the interpolating situation of partial coupling, some approxi- 
mate analytical results have been derived. In Ref. [2l[, for instance, it has been found 
that the density p(E) for arbitrary A and m is given by 

m M 

8=1 

where 

' 2 -v^TTP, \E\<a 
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is Wigner's semi-circle law with a = ^jN/ot and 
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The transition parameter utilized in the following is defined as 22| 



(9) 



A = A 2 p(0) 2 = A 2 -^-. (10) 

Eq. (jSJ) can be used to calculate exactly analytically the NND for 2x2 and 3x3 
matrices [6j. For the 2x2 case the DGOE, Eq. (jSJ), gives 

P 2 , 2 ( S ,A) = ^/o(^)exp(-^), (11) 



where Iq is the modified Bessel function, whose asymptotic form is 

I (x) - -^=. (12) 

Thus, ^2x2(5, 0) = 1, and there is no level repulsion for A — > 0. In the opposite limit, 
A — > 00, I (x) ~ 1 — x 2 /4 and one obtains 

P 2x2 ( S ,A»l)^ S ^ e - 2 /16A. (13) 

For higher dimensions Eq. (jSJ) can only be used for numerical simulations, however, using 



appropriate perturbative methods Leitner [23] was able to find a formula for the NND. 



He started basically with the formula for the nearest neigbhour spacing distribution for 
the superposition of m GOE's block matrices [l[ 

d 2 

where, for the case of all block marices having the same dimension one has 

E m (s)= f^i(-)) m , (15) 



m 



POO 

E x [x)= / (1-F(t))dt, (16) 

J X 



F(t)= / P l (z)dz. (17) 





In the above Pi(z) is the normalized nearest neighbour spacing distribution of one block 
matrix. It is easy to find for P m (s), the following 

P m (s) = - ^{s/m))^ 1 P 1 {s/m) + {m-l){E 1 {s/m)r-\l - F{s/m)f] (18) 
m 

= PgX8)+P%>(8) (19) 

If all the block matrices belong to the GOE, then one can use the Wigner form for P\{z) 

P 1 (z) = \ze-^ « *-z, (20) 



thus 

F1 (z) = \-er-^ w 

4 



F 1 (z) = l-e-^ z2 w -z 2 , (21) 



E 1 (z) = erfc\^z) « 1 - z. (22) 

where the large- z limits of Eqs. (I20p - (I22Q are also indicated above. It is now clear that the 

above expression for P m (s), ( |T8l) and ( |T9l) . contains a term Pm\s) with level repulsion, 
indicating short-range correlation among levels pertaining to the same block matrix and 
a second term Pm\s) with no level repulsion, implying short-range correlation among 



NND levels pertaining to different blocks. Notice that for very small spacing, P m (s) 
behaves as 

P m (s)*^- 2 s + — (23) 

for m = 1, we get the usual -Pi(O) = 0, while for m > 1, we get P m (0) — (m — l)/m. 
To account for symmetry breaking, Leitner 23j considered the mixing between levels 
pertaining to nearest neigbhour block matrices. This amount to constrain the mixing to 
be of the form given by Eq. (11 II) for the 2x2 DGOE, thus, he found 



P m (s, A) = P«( a ) + P 2x2 (s, A)P^(s). 



(24) 



Though P m (s) is normalized, P m (s, A) is not. Accordingly one supplies coefficients 
and Cd, such that 

P m (s, A, c^, c d ) = c N P m (c D s, A) (25) 

is normalized to unity. Similarly, < s > should be unity too. Eq. fl24|) can certainly 
be generalized to consider the effect of mixing of levels pertaining to next to nearest 
neighbour blocks, and accordingly, P 3x3 (s,A), given in Ref. |6j would be used in Eq. 
(124)) instead of P2 X 2(s, A). In the following, however, we use Eqs. (|24|) . fl25|) as Leitner 
did 23 . 



In Ref. [23j, Leitner also obtained approximate expression for the spectral rigidity 
A 3 (L) using results derived by French et al. (24J. Leitner 's approximation to A 3 is equal 
to the GOE spectral rigidity plus perturbative terms, that is 



A 3 m) (L;A) 



A 3 (L;oo) + 



m — 1 



1 



e 2 L 2 



x ln(l + e 2 L 2 ) + — tan _1 (eL) 
eL 



2e 4 L 4 
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where 



TT 



2(r + tt 2 A) 



For the cut off parameter we use the value [25j r 



Cm & 



7 1 , where c r , 



(26) 



(27) 



m m/(m-l) 



and 7 ~ 0.5772 is Euler's constant. This choice guarantees that when the symmetry is 
not broken, A = 0, A^ m) (L,0) = mA 3 (L/m, oo). In Ref. 0, Leitner fitted Eq. (ESD 



for m=2 to the NND from Ref. [17j, however, he did not fit the spectral rigidity. It is 
often the case that there are some missing levels in the statistical sample analysed. Such 
a situation was addressed recently by Bohigas and Pato [27| who have shown that if g 
fraction of the levels or eigenfrequencies is missing, the A 3 (L) becomes 



Af(P) 



9 



L 
15 



S) 2 A 3 



L 



1-9 



:28) 



The presence of the linear term, even if small, could explain the large L behavior of the 
measured A 3 (L). We call this effect the Missing Level (ML) effect. Another possible 
deviation of A 3 from Eq. fT2"6"|) could arise from the presence of pseudo-integrable effect 



(PI) |20L |28I |. This also modifies A 3 by adding a Poisson term just like Eq. 



We now apply our model to analyse the eigenfrequency data of the elastomechanical 
vibrations of an anisotropic quartz block used in [13]. In this reference in order to break 
the flip symmetry of the crystal block gradually they removed an octant of a sphere of 
varying size at one of the corners. The rectangular quartz block has the dimensions 14 x 
25 x 4 mm 3 . The radii of the spheres containing the octants are r = 0.0, 0.5, 0.8, 1.1, 1.4 
and 1.7mm representing figures (a) — (/). Figs, lx and 2x of Ref. l3| correspond 



to the octant r > 1.7. They found 1424, 1414, 1424, 1414, 1424 and 1419 frequency 
eigenmodes, respectively. The histograms and circles in the two figures of Ref. [17| 
represent the short-range nearest-neighbor distributions P{s) (Fig. 1) and the long range 
A 3 (L) statistics (Fig. 2). 

The results of our analysis are shown in the two figures. In Fig. 1, the sequence of 
six measured NNDs were fitted for m = 2 and m = 3. It can be seen that the DGOE 
model with three coupled GOE's give a comparable and in some cases even better fit 
than the m = 2 one. Figure la in fact shows a rather sharp peak in our calculated 
P(s) for m = 3, -P^s, 0.0056). We consider this a failure of our formula ([25]) for the 
uncut crystal. In fact, a more appropriate description of the uncut crystal is to take 
A = 0, namely a superposition of 3 uncoupled GOE's, which works almost as good as 
the 2 uncoupled GOE's description. The other parts of figure 1, (o) — (x) seem to show 
the same insensitivity of P m (s,A) to m; the number of matrix blocks used in DGOE 
description. It is this insensitivity of the short-range nearest neighbour level correlation, 
measured by the spacing distribution, to the assumed symmetry inherent in the uncut 
crystal (and thus the number uncoupled GOE's employed to describe it) that forces us to 
examine the long-range level correlation, namely spectral rigidity, "measured" by Dyson's 
A 3 statistics. 

In Fig. 2. the A-statistic was fitted with equation f[2"6"j) . It is clear from the figure that 
a good fit to the data of Ref. [17| is obtained with m = 3 for the values of A given in table 
1. This is to be contrasted with the case of m = 2 which, according to Eq. f[2"6l) results 
in Af\L, A) that is always below the one with A = 0, A3^(L,0), which itself is always 
below the data points of Ref. [13]. For this reason, only the A 3 2 ^(L,0) is shown in the 
figure. It should be noted that the A-statistics of the uncut crystal, Fig. 2a is very well 
described by that of 3 uncoupled GOE's, namely Ag (L) = 3A 3 1 ' ) (L/3) which is always 
larger than the above mentioned A^(L) = 2A 3 1 ^(L/2). The most conspicuous exception 
is Fig. [2b which corresponds to r = 0.5 mm and where 1414 frequency eigenvalues were 
found. We consider this a potential ML case and take for A3, the expression given in 
Eq. (1281) and use it in Eq. ([26]). We find perfect fit to the data, if g is taken to be 6%, 
namely only 94% of the eigenfrequencies were in fact taken into account in the statistical 
analysis. In contrast, if 2GOE is used we still do not get very good agreement even if 
18% of the levels are taken to be missing, as shown in Fig [3J There is, threfore, room 
to account much better for all cases (Fig. 2a, 2c, . . . ) by appropiately choosing the 
correponding value of g. 

In conclusion, a random matrix model to describe the coupling of m-fold symmetry 
is constructed. The particular threefold case is used to analyse data on eigenfrequencies 
of elastomechanical vibration of a anisotropic quartz block. By properly taking into ac- 
count the ML effect we have shown that the quartz block could very well be described 
by 3 uncoupled GOE's , which are gradually coupled by the breaking of the three-fold 
symmetry (through the gradual removal of octants of increasing sizes), till a 1GOE situa- 
tion is attained. This, therefore, indicates that the unperturbed quartz block may posses 



TABLE I: Values of A obtained by fitting Eqs. (I25D and (I26h respectively to the experimental 



NNDs and spectral rigidities from Ref. 171 ]. 





POO 


A 3 (L) 


Data Set 


Ref. [26] 


Eq. (j2U m=2 


Eq. {25]) m=3 


Eq. (|25|| m=3 


(a) 


0.0013 


0.0030 


0.0067 


0.0056 


(b) 


0.0054 


0.0063 


0.0098 


0.0016 


(c) 


0.0096 


0.010 


0.017 


0.0017 


(d) 


0.0313 


0.032 


0.064 


0.027 


(e) 


0.0720 


0.070 


0.13 


0.050 


(f) 


0.113 


0.12 


0.30 


0.16 


(x) 


0.138 


0.13 


0.34 


2.4 



another symmetry, besides the flip one. We have also verified that if a 2GOE description 
is used, namely, m — 2 , then an account of the large-L behaviour of A3 can also be 
obtained if a much larger number of levels were missing in the sample. In our particular 
case of Fig. 2b, we obtained g = 0.18. This is 3 times larger than the ML needed in the 
3GOE description. We consider the large value of g needed in the 2GOE description, 
much too large to conform to the reported data in Ref 17| . A preliminary version of the 
formal aspect of this work has appeared in j29|. 
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FIG. 1: Nearest Neighbour Distributions. Histograms show data (a)-(x) from Ref. [13]. Thin 
and thick solid lines show fits to the data carried out using Eq. (|25[) with m=2 and m=3 
respectively. In graph (x) the long-dashed line is the Poisson distribution, the dot-dashed line 
is the Wigner distribution and the dashed and dotted lines are the respective distributions for 
superpositions of 2 and 3 uncoupled GOEs. See Table I for the values of A obtained from the 
fits and the text for details. 




FIG. 2: Spectral Rigidities. Circles show data (a)-(x) from Ref. [17]]. Thick solid lines show 
fits to the data carried out using Eq. ([25]) with m=3. The dot-dashed line is the GOE spectral 
rigidity and the dashed and dotted lines are the respective rigidities for superpositions of 2 and 
3 uncoupled GOEs. See Table I for the values of A obtained from the fits. 




FIG. 3: The Missing Level Effect: The full curve is the 2G0E A 3 for case 2b with 18% of the 
levels missing, while the dashed curve corresponds to 3G0E with only 6% of the levels missing. 
The data points (open circles) represent fig. 2b. See text for details. 



